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A new formulation of piezo-resonator theory, called for 
convenience the ‘‘mechanical theory,”’ is developed for a 
crystal resonator with gap between crystal and electrodes, 
involving the increase in effective stiffness caused by space 
or surface charge, and explaining the increase in resonance 
frequency with increasing gap. The cases treated are 
thickness vibrations (both compressional and shear) of a 
plate, and lengthwise compressional vibrations of a bar. 
Specialized values of the elastic, dielectric and piezoelectric 
constants are employed, depending upon the nature of the 
elastic reactions for each cut and for each mode of vibra- 
tion. The mechanical theory is compared with the cus- 
tomary ‘‘electrical theory,”’ according to which the resonator 
maintains a constant stiffness, the increase in frequency 
with increasing gap being accounted for by the electrical 
reactance of the gap considered as a series capacitance. 
When the crystal surfaces are equipotential (thickness 
vibrations, or lengthwise vibrations of a bar whose surfaces 
have been covered with a metallic deposit), the two 
theories are in agreement. It is shown that the mechanical 
theory, unlike the electrical theory, describes correctly 
the performance of a vibrating crystal bar that is not 
metallically coated. The characteristics of the resonator 
and the contrast between the two theories are illustrated 
by resonance-circle and phase diagrams. The relative 
increase in frequency caused by the gap is shown to be 
Af/fo= Uw/(e+kw), where fo is the frequency at zero gap, 


fot+Af the frequency with a gap of width w, e the thickness 
of the crystal in the direction of the applied electric field, 
k the dielectric constant and U=16e/mc, « being the 
appropriate piezoelectric constant and c the effective 
elastic constant at zero gap. In the case of the uncoated 
bar, the function U has to be multiplied by ?/8. Observa- 
tions were made with eight quartz crystals of the following 
cuts: Y-cut, X-cut, BT- and AC-cuts (notation of Lack, 
Willard and Fair) for thickness vibrations, and an X-cut 
bar for lengthwise vibrations. With small gaps, the experi- 
mental values of the function U are, except with the BT- 
and AC-cuts, from 13 to 60 percent smaller than the 
theoretical values. In most cases the discrepancy becomes 
still greater at large gap values. The probable reasons for 
these discrepancies are discussed, with the conclusion that 
a large increase in frequency with gap may be taken as a 
criterion of a good resonator (large piezoelectric constant, 
freedom from twinning and coupling). A Rochelle salt 
Y-cut 45° bar when tested for the gap effect showed a 
satisfactory agreement with theory, considering the 
uncertainty in the values of the dielectric and piezoelectric 
constants. It is shown that the effect of gap upon frequency 
can be utilized in determining approximately the piezo- 
electric constant of a crystal. The theory of the dependence 
of frequency upen gap for overtone vibrations is briefly 
discussed. 


INTRODUCTION 


IEZOELECTRIC resonators are frequently 

provided with a free space between the 
crystal and one or both electrodes, although 
there is a growing practice of using as electrodes 
thin metallic films deposited directly on the 
crystal surfaces. When present, the gaps modify 
the over-all mechanical and electrical constants 
of the resonator ; they thus affect both amplitude 
and natural frequency of vibration. The gap 
effect is dependent on space and surface charge 


phenomena associated with the cut and mode of 
vibration employed. It yields to simple theo- 
retical treatment only in the cases of relatively 
thin plates and of relatively long rods. Even with 
this simplification, a rigorous theory of the 
resonator, with or without a gap, would have to 
take account of boundary conditions, all elastic, 
frictional and electric constants of the material, 
coupling effects between vibrations of various 
modes, nonuniformity of applied electric field 
and, when air gaps are present, resonance effects 
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in the air itself. Solutions of special problems 
involving certain of these considerations have 
been offered by many writers. 

The present paper is restricted to thickness 
vibrations of plates and lengthwise vibrations of 
bars, modes which are of both scientific and 
engineering interest. The present treatment 
differs from previous discussions chiefly in offer- 
ing a more complete explanation of the influence 
of the gap upon the effective stiffness of the 
crystal. It is assumed that only one degree of 
freedom need be considered ; that the impressed 
electric field remains uniform even with the 
largest gaps; and that frictional and resonance 
effects due to air in the gap may be ignored. 
The discrepancies which these idealized con- 
ditions introduce between theory and observa- 
tion will be discussed in §19. 


I. THEORY OF THE RESONATOR 
The electrical theory 


$1. It has become customary to consider the 
crystal as equivalent to a pure dielectric capaci- 
tance C,, Fig. 1a, in parallel with the series chain 
RLC.'.? When a gap is present, it is treated as 
a series capacitance C2.?:* R, L, C and C, are 
regarded as fixed constants for each particular 
resonator, as are also the mechanicalconstants 


R ¢ 
Qa 
i, 
R' 


Fic. 1. Equivalent electrical constants of crysta! with gap. 


'K. S. Van Dyke, Phys. Rev. 25, 895 (1925); Proc. 
I. R. E. 16, 742 (1928). 

2D. W. Dye, Proc. Phys. Soc. London 38, 399 (1926). 

%Y. Watanabe, Elec. Nach. Tech. 5, 45 (1928); Proc. 
I. R. E. 18, 695 (1930). 
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of the crystal. We shall for convenience refer to 
this theory as the “electrical theory,” in contrast 
to the ‘‘mechanical theory’’ which is about to be 
developed. Watanabe* has shown that for the 
network in Fig. la can be substituted, at all 
frequencies, that shown in Fig. 1b, where 


(1) 
Cc’ 2 C2 
Ciy? C Ci 
- (1+) +-(14+2), (2) 
Ce Cy C2 
=1+Ci/C2. (3) 


R, L, C and C, are the equivalent electrical con- 
stants for the crystal with zero gap, while C; is 
the capacitance of the gap. Equivalent networks 
of substantially the same type had previously 
been treated by Johnson and Shea.‘ Although 
the equivalent electrical constants have fre- 
quently been measured, few writers have paid 


any attention to the effect of the piezoelectric: 


polarization field on the elastic constant, and 
none, so far as is known to the writer, has con- 
sidered the effect of the gap on this quantity. 
It has been customary to treat the elastic con- 
stant as a fixed quantity, unaffected by the gap. 

It is already known that space and surface 
charges caused by polarization of the crystal 
may have effects on the resonator constants,>- 
but hitherto this consideration has been applied 
only to a few special cases, in none of which has 
the universal validity of Eqs. (1) and (2) been 
questioned. In the present paper these effects are 
treated more rigorously and with greater gener- 
ality, in order to see under what circumstances 
Eqs. (1) and (2) hold. Thickness vibrations of 
thin plates, of compressional and shear modes, as 
well as longitudinal vibrations of bars, are con- 


4 5} S. Johnson and T. E. Shea, Bell Sys. Tech. J. 4, 52 
(1925). 

5 E. Giebe and A. Scheibe, Ann. d. Physik 9, 93 (1931). 

®° 1. Koga, Physics 3, 70 (1932). 

7 J. Kobzarew, J. Appl. Phys. (Moscow) 6, 17 (1929). 

5M. v. Laue, Zeits. f. Physik 34, 347 (1925). 

aa) Giebe and A. Scheibe, Zeits. f. Physik 46, 607 
(1928). 

10 D) Koga, Rep. Rad. Res. and Works in Japan 2, 157 
(1932). 
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sidered. Application is made to several well- 
known cuts in quartz, and the results are com- 
pared with experimental data. 


The mechanical theory 


(a) Thickness vibrations 

§2. The generalized theory for thin plates will be 
developed first, applicable to all crystals and all 
cuts which satisfy the requirement that an 
alternating electric field normal to the surface of 
the plate causes vibrations (compressional or 
shear) with wave propagation parallel to the 
field. For convenience, this normal direction will 
here be denoted by x, and elastic strains, whether 
of compressional or shear type, by x,; when 
application is made to particular cuts, the axial 
direction appropriate to each cut will be sub- 
stituted (Table I). Strictly, the theory is appli- 
cable only to plates of infinite area, in which case 
the only strain is that characteristic of the thick- 
ness vibration. If the boundaries are finite, then 
as is well known various other vibration modes 
may become excited, due not only to primary 
piezoelectric effects but to various elastic re- 
actions (coupling effects) as well. These other 
modes may be expected to complicate the 
electric conditions in the crystal, and thus to 
modify the conclusions derived from the present 
theory, as will be discussed in §19. It would be 
extremely difficult to include these disturbing 
effects in the theory, and perhaps even more 
dificult to test so complete a theory experi- 
mentally. In any case, even if coupling effects 
were taken into account, failure of experimental 
results to conform very exactly to theory must 
be expected, owing, even when great care is 
expended, to such unavoidable causes as non- 
uniformity of electric field, nonparallelism of 
crystal and electrode surfaces, imperfect orienta- 
tion of crystal plate, mechanical defects in the 
crystal, and especially, when adjustable elec- 
trodes are used, the impossibility of taking 
readings with the gap at strictly zero value. 

The resonator is a plate of piezoelectric 
crystal of thickness e, between conducting elec- 
trodes, as shown in Fig. 2. The surface area is 
assumed to be so great that boundary effects can 
be ignored. The total gap is w; under our simpli- 
fying assumptions it does not matter whether 
the gap is all on the same side of the crystal or 
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Fic. 2. Resonator for thickness vibrations. 
\ 

not. The possible effect of the position of the 
crystal between the electrodes on the frequency, 
for a finite plate, is discussed in §20. The in- 
stantaneous driving voltage impressed on the 
electrodes is V= Vo cos wt, and hence the driving 
voltage across the crystal is 


Vi=eV/e’, (4) 


where e’ =e+kw. For k we take the value of the 
dielectric constant appropriate to the orientation 
of the particular crystal cut used. (As a matter 
of dimensional consistency it should be said that 
k is here strictly the ratio of the dielectric con- 
stant of the crystal to that of the gap.) V; causes 
in the crystal an instantaneous driving field 
E,=V,/e=V/e’. The corresponding field in the 
gap is 

(5) 


Since we are dealing with a one-dimensional 
problem and are -confining ourselves to the 
resonance region for the fundamental thickness 
frequency, we may substitute for the crystal an 
equivalent mechanical system of concentrated 
mass M, frictional resistance W, and stiffness 
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including effect of gap G’. This method was first 
used for lengthwise vibrations": and has now 
become common practice. The equation of 
motion, for the resonance region, is 


ae dé 
M—+W—+G’'t= Fy cos wt, (6) 
dt? dt 


where & is the equivalent mechanical displace- 
ment, equal to the actual displacement of either 
face of the crystal, the nodal plane being at rest. 

A word should be said at this point concerning 
shear vibrations in y-cut quartz plates. The same 
considerations apply equally to all types of 
vibration in which the axis of shear lies in the 
plane of the plate, and this includes the AC and 
BT cuts treated in §8. For the y-cut quartz 
resonator the stress and strain are X, and x,, 


respectively, and the direction of & lies in the 


plane of the plate, instead of being normal to 
the surface as in the X-cut. In the neighborhood 
of resonance the distribution of strain is sinu- 
soidal, exactly as with compressional waves. 
The equivalent mass M is the same, as may be 
proved by consideration of kinetic energy. X in 
Eq. (8) means in this case a tangential stress X ,. 
All the resonator equations here given, and 
hence the equivalent electrical networks, are 
applicable equally to shear and compressional 
modes. Numerical values vary, of course, from 
one type to another, as shown in Table I. 

As in the case of a longitudinally vibrating 
rod," the equivalent instantaneous sinusoidal 
impressed force (positive when extensional), is 


F=F,y cos w= —2AX =2AE; 
V/e’=(2eA Vo/e’) cos wt, (7) 


where X is the instantaneous driving mechanical 
stress for thickness vibrations, given by the 
equation 


X= (8) 


in which ¢ is the piezoelectric constant relating 
Ey to X. 

According to Voigt’s terminology, which is 
here followed, a stress is positive when com- 
pressional, while a strain is positive when 


"W. G. Cady, Proc. I. R. E. 10, 83 (1922). 
2 W.G. Cady, Phys. Rev. 19, 1 (1922). 


extensional. This, together with the assumption 
that the driving force F is positive when com- 
pressional, accounts for the negative signs jn 
Eqs. (7) and (8). 

On solving (6) and substituting from (7) we 
have for the maximum displacement (using 
primed quantities for the general case where a 
gap is present) 


Fo 2AVy 22AVo 
= = Y,’ (9) 
wZ, we'Z,! we’ 


The mechanical impedance of the crystal js 
given by 
Z.° =W?+(wM —G’'/w)?. (10) 


Y.’=1/Z.’ represents the mechanical admit- 
tance. The displacement at time ¢ is 


sin (wt—@), (11) 


where —ttan 0=(wM—G'/w)/W. (12) 


For the velocity of motion at the surface of 
the crystal we find from Eq. (11) 


v=dt/dt=wk cos (wt—6). (13) 


From Eq. (9) it is clear that v is a maximum at 
that frequency for which Z, is a minimum, the 
mechanical reactance being then zero. 

The equivalent mechanical constants, for a 
plate of large area A, are found, by the same 
method as for lengthwise vibrations,’ to be 


M=peA/2, W=7°pQA/2e, G’=2°c'A/2e, (14) 


where p and Q are the crystal density and vis- 
cosity, and c’ the effective elastic constant 
derived in §4. As will be seen, G’ and c’ them- 
selves are functions of the gap. The equation for 
G’ is derived from that for M together with the 
resonance relation Mw=G’/w and the fact that 
for the fundamental vibration frequency (c’/»)! 
=2fe. In actual resonators, W usually depends 
chiefly on the mounting and surrounding air. 
The true quartz viscosity Q is very small; no 
reliable values for thickness vibrations seem to 
have been published. Van Dyke" finds for an 
X-cut Y-wave quartz bar vibrating lengthwise 
in vacuum a value of the ‘internal friction 


8K. S. Van Dyke, Proc. I. R. E. 23, 386 (1935). 
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constant” (ratio of logarithmic decrement to 
Young’s modulus) equal to 7 (10-'), whence 
from the equations in §6 the value of Q is 1.1. 
If the decrement of the resonator is measured, 
the over-all value of W may be found from the 
equation A=W/2fM, f being the resonance 
frequency and A the logarithmic decrement. 

In order to deduce the equivalent electrical 
constants, we seek an equation for the in- 
stantaneous current J flowing to the resonator. 
This is of two components (Fig. 1): 
(15) 
C;' represents the pure dielectric, or ‘‘geo- 
metrical,’ capacitance of crystal and gap to- 
gether, including the stray capacitance of 
mounting and leads. It is an important param- 
eter and should preferably be measured directly 
for each gap whenever the theory is applied to 
experimental results. Theoretically, ignoring 
stray capacitance and assuming the crystal 
surfaces to be equipotential, C,’ is the capaci- 
tance of two condensers C,=kA/4ze and C2 
=A/4rw in series (Fig. 1), so that C;’=kA/4ze’ 
=C1C2/(Ci+ C2), as in Eq. (3). C, had best be 
measured rather than calculated. Unless the gap 
is relatively large, the above equation for C2 can 
safely be assumed to be approximately correct, 
and this assumption is made throughout (see 
Sec. 19).14 

§3. Equivalent electrical constants. The second 
component of current, i, in Eq. (15), which may 
be called the piezo-current, owes its existence to 
the piezoelectric displacement D. This quantity 
is derived as follows. We assume that the strain 
distribution in the crystal is sinusoidal in both 
space and time, and that at any given time / 
the equation for the instantaneous strain at 
distance x from the central, or nodal, plane of 
the crystal (for frequencies in the neighborhood 


_“In the author's original paper, reference 11, that por- 
tion of the deformation of the crystal that would have 
been produced by a static potential difference V was re- 
ferred to as the “equilibrium elongation.” This complica- 
tion is now aveiiel by taking for the dielectric constant 
a value based upon the constraints that characterize the 
particular type of vibration (§10). Coupled modes are here 
ignored. A further departure from the original paper lies 
in defining the electric displacement in the dielectric as 


=kE, in conformity with modern usage, instead of 
D=kE/4r. 
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Fic. 3. Distribution of polarization P, space charge p, 
electric intensity E and potential V, in a resonator for 
thickness vibrations. Crystal thickness AB =e. 


of the fundamental) is 
X2=Xz29 Cos (rx/e) 


where xz is the strain at the center. The origin, 
as shown in Fig. 3, is at the center of the plate 
AB. The gap is w=A’A+BB’. The curve for x, 
is not shown, but it would be similar to the 
curve P. 

In order to establish a relationship with Eq. (9) 
and thence to find the equivalent electrical 
constants, it is desirable to express D in terms 
of &, the instantaneous mechanical displacement 
at point B in Fig. 2, at which x=e/2. It is easily 
verified that x,9=7£/e, whence 


cos (rx/e). (16) 
At any x the piezoelectric polarization is 


TE TX TX 
P(x) =ex,=—cos—=Pycos—, (17) 
e e e 


where Po=zet/e. This gives rise to a space 
charge 


dP(x) wPo 
p(x) = ————_ =——- sin —=posin—, (18) 
dx e 


where 


The space charge causes a component of 
electric intensity E shown in Fig. 3, the value 
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of which at any x in the crystal has been derived 
in a former paper :'° 


4poe 


(19) 


The second term in Eq. (19) is that which 
results from writing E(x) = —4aP(x)/k directly. 
Koga*:' derives an expression equivalent to 
this, but he does not consider the first term, 
which depends on the distribution of P and 
expresses the gap effect. The curve V in Fig. 3 
represents the potential distribution in crystal 
and gaps, on the assumption that the electrodes 
are at zero potential. The potential equation is 
given in reference 15, but it is not needed in the 
present paper. 

The total field at any point x in the crystal 
has two components, E(x) due to the strain, 
and the impressed field E,= V/e’ from Eq. (1). 
By deriving E(x) on the assumption of zero 
potential difference between electrodes we have 
separated the two components. FE, accounts for 
the first term in Eq. (15) and has nothing to do 
with the piezoelectric reactions. The second 
term in Eq. (15) is derived from E(x), as will 
now be proved. 

By setting x= +e/2 in Eq. (19) we find for the 
electric intensity in the gap due to the strain, 
and therefore for the corresponding displacement 
component, 

Ey = D=8ret/e’. (20) 


It may be noted that the displacement D, 
which by the principle of continuity is the same 
throughout the resonator, depends entirely on 
the value of E at the crystal surfaces. The effect 
of the form of the E curve on the mechanical 
stiffness will be discussed in §4. 

From Eqs. (7), (15) and (20) the resonator 
current is 


kA dV dg 
4re’ dt e’ dt 
kwA V9 42A?V 
sin wt-+— cos (wt—@). (21) 
4re’ eZ. 


W. G. Cady, Physics 6, 10 (1935). 
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We are concerned at present only with the 
second term, which represents the piezo-current 
7 flowing through the equivalent electrical im- 
pedance denoted by Z’. The maximum current 
value can thus be expressed as ip= V)/Z’ 
= Vo/|R"°+(wL’—1/wC’)*}!, where R’, L’ and 
C’ are the equivalent electrical constants. Com- 
parison of this equation with Eqs. (10) and (21) 
shows at once that, calling 


=e" /42A?, (22) 
Z'=r'Z,, R’'=r'W, L’=r'M, C’=1/r'C’. (23) 


The primes refer to crystal with gap. When there 
is no gap, r’ reduces to r=e?/4e°A? and Eq. (23) 
becomes 


Z=rZ., R=rW, L=rM, C=1/rG. (24) 


These values of R, L and C are in agreement with 
those first derived by Van Dyke.' Just as the 
crystal without gap has as its equivalent elec- 
trical constants the capacitance C; in parallel 
with the series chain R, L, C, so when a gap is 
present, the equivalent constants consist of C/’ 
in parallel with R’, L’, C’. This is evident from 
Eq. (15) or (21), the two terms of which represent 
the current in the parallel condenser C,’ and the 
piezo-current in the R’L’C’ chain, respectively. 

Comparison of Eq. (23) with (24) shows that 
the effect of the gap upon the equivalent elec- 
trical constants is contained in the electro- 
mechanical ratio 7’, which involves the quantity 
e’=e+kw, w being the gap width ; and, moreover, 
in the case of C’, the mechanical stiffness G’, 
when the gap is present, is slightly greater than 
G for zero gap, as will presently be explained. 
Consideration of the ratio r’ shows that the 
combination of crystal and gap is electrically 
equivalent to a crystal of the same dielectric 
constant with adherent electrodes, of thickness 
e’, having its mechanical resistance and mass 
increased in the ratio (e’/e)?, and its stiffness in- 
creased in a slightly higher ratio. 

$4. Equations for elastic constants. We turn now 
to the effect of the gap upon the equivalent stiff- 
ness G’ and effective elastic constant c’, the 
relation between which was expressed in Eq. 
(14). Consider the reacting stresses when the 
crystal surfaces are displaced by an amount { 
at e/2, Fig. 3 and —£ at —e/2 (points B and A, 
respectively). These reacting stresses do not 
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include the driving stress X of Eq. (7). For a 
given £, we have seen that the strain distribution 
js given by Eq. (16): 


x, =(mt/e) cos (xx/e)=X,/c’. (25) 


In the last member of Eq. (25) the strain at 
any point is expressed in terms of total stress 
(exclusive of the driving stress) and effective 
elastic constant. The chief component of X;, is, 
of course, the pure elastic stress, which for 
generality we will call Xo=—coxz, co being the 
pure elastic constant appropriate to the material 
and orientation of the plate. For example, with 
an x-cut quartz plate this takes the familiar 
fom X,=—cyx,. In a non-piezoelectric sub- 
stance Xo is the only reacting stress. When the 
material is piezoelectric and under uniform 
strain in a direction parallel to the field (length- 
wise vibrations, §6), then with zero gap Xo is 
still the only reacting stress; when a gap is 
present, there is an additional stress component 
due te polarization surface charge. With thick- 
ness vibrations the strain is sinusoidally dis- 
tributed, and the resulting space charge, as will 
be shown, causes two additional stress com- 
ponents X, and X¢, the latter of which is depend- 
ent on the gap. We now proceed to derive the 
elastic constants, or, as they may more appro- 
priately be called, the quasi-elastic constants, 
corresponding to these two stress components. 

From Eq. (19) and Fig. 3 it is seen that the 
electric intensity E due to space charge is of 
sinusoidal form, having negative values in the 
interior of the crystal and positive near the outer 
surfaces. The vertical displacement of the curve 
is expressed by the first term in Eq. (19), which, 
by the factor e’, depends upon the gap. The sig- 
nificance of this becomes evident when we con- 
sider the mechanical stresses that E causes in the 
crystal, through the converse piezoelectric effect, 
which in generalized form may be written X =cE. 
On applying this equation to Eq. (19) there 
results 


mx 2e 
cos —), (26) 
e 


where 


and X2=87et/ke’. (27) 
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It is easily verified that wherever E is positive, 

the resulting stress has the same sign as the 

assumed strain x,; it therefore is opposed to the 

pure elastic restoring stress and tends to lower 

the vibration frequency. Where E is negative, as 

is predominantly the case, the opposite is true. 
From Eqs. (16), (26) and (27), 


— (402 /k)x,= (28) 


The stress X, at any point is thus seen to be a 
constant times the strain at the same point, 
whence it is evident that c, the coefficient of 
xz, is of the nature of an elastic coefficient that 
is constant throughout the crystal. This quasi- 
elastic constant has also been derived for certain 
special cases by other writers.°: © Together with 
the constant ¢2, it should never be overlooked 
when elastic constants are computed from 
thickness vibrations. For example, it offers an 
explanation of the discrepancy between theory 
and observation noted by Bechmann.'® Accord- 
ing to Voigt’s notation, the negative sign in Eq. 
(28) indicates that c, is to be added to the pure 
elastic constant Co. 

The case is not so simple with Xe, as it is not 
a sinusoidal term. It is rather a stress uniform 
throughout the crystal and is therefore to be 
treated in a manner analogous to that used for 
the driving stress X, Eq. (7). Like X1, Xo is pro- 
portional to and in phase with &, and hence is not 
a part of the driving stress. Nevertheless, it may, 
like X in Eq. (7), be regarded as a stress applied 
at the surfaces of the resonator. We may therefore 
regard it, from the point of view of the equivalent 
mechanical system with concentrated mass M, 
as supplying a component of the total restoring 
force G’ acting on M when the instantaneous dis- 
placement is &. As was shown for the case of X 
in previous papers," we must take 2AX_ as 
this component. We thus obtain the equivalent 
stiffness component G2 from Eq. (27): 


= 
where — /ke’. (29) 


The negative sign results from Voigt’s conven- 
tion respecting the definitions of stress and 
strain. 


16 R. Bechmann, Zeits. f. tech. Physik 16, 525 (1935). 
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To G, there corresponds the second of the 
quasi-elastic constants, ¢2. The vibration fre- 
quency is determined by the three quantities Co, 
c; and ¢2, whose sum is the total effective elastic 
constant c’. Their relation to the equivalent 
stiffness G’ may be expressed by writing Eq. 
(14) thus: 


G’ = /2e) (cotcit+ce). (30) 


The chief component of stiffness is Go, related 
to the pure elastic constant co by the equation 
Go=7’°Aco/2e. From (30) it is seen that G, 
= 7*Ace/2e, whence from Eq. (29) 


—32e€e/rke’. (31) 


From Eqs. (30) and (31), together with (28), the 
final form of the effective stiffness G’ is found 
to be 


G’ =Go+Git+G, 


4re 32€e rAc’ 
2e 2e 
At zero gap G’ becomes 
TA 32é wAc 
G=— = (33) 
2e k wk 2e 


The effective elastic constant c’ for a thickness- 
vibration resonator with gap is therefore 


(34) 


This equation shows that the effective elastic 
constant c’ is at zero gap slightly greater than 
the pure elastic constant co and that the value 
increases with increasing gap. Herein lies the 
explanation of the change in the stiffness G’ with 
gap referred to above. The gap effect is contained 
in the third term, which approaches zero for 
infinite gap. 

The values of the piezo-capacitances C’ and C 
in terms of the fundamental constants, from 
Eqs. (22), (23), (24), (32) and (33) are C’ 
and C=8eA/n’ec. 

For an X-cut quartz plate vibrating in the 
thickness mode the correction term c, has a 
value about 0.87 percent as large as co, while for 
zero gap C2 is about 0.71 percent as large as Co. 
Thus for zero gap the effective elastic constant 


c’, for this cut, through the combined effects of 
c, and ¢2, becomes about 0.16 percent greater 
than the pure elastic constant co. The natural 
resonator frequency, as the gap increases from 
zero to infinity, would therefore be expected to 
rise by about 0.3 percent (see §19). 

The fact that G’ in Eq. (32) consists of three 
terms suggests that C’ in Eq. (23) may be re- 
garded as equivalent to three capacitances in 
series. The values of the three are easily derived, 
but are of no importance in the present dis- 
cussion. 

We now consider briefly the effect of varying 
gap upon the curves in Fig. 3. For the same 
amount of strain, the curves for P and p remain 
unaltered, but the E curve, by Eq. (19), moves 
vertically downward as the gap increases, while 
the abscissas for maximum and minimum YV 
shift correspondingly. Let the values of x cor- 
responding to E=0 (maximum and minimum /) 
be denoted by +x». It is easily shown that, for 
gap W, COS TX,/e=2e/re’ where e’=e+kw as 
before. When w=0, cos tXm/e=2/2, whence 
Xm =0.556e/2, so that the E curve cuts the axis 
at points nearly midway between the center and 
the surfaces of the crystal. When w=, 
COS TXm/e=0 and x,=+e/2. The E curve then 
touches the points A and B without intersecting 
the axis at all, which corresponds to the vanishing 
of the third term in Eq. (34) at infinite gap. 

§5. Equivalence of the electrical and mechanical 
theories. In the resonator theory given in §§2 to 4, 
the crystal with gap was treated as a unit, with 
a mechanical stiffness G’ dependent on the gap. 
A necessary postulate, which was also implicit 
in the “‘electrical’’ theory, is that each surface of 
the crystal be equipotential; this condition is 
satisfied so long as the only wave-propagation is 
in the thickness direction. 

It has already been stated that the values of 
R, L and C.in Eq. (24), derived by the method 
of §3, are in agreement with Van Dyke’s values, 
which are also used by Watanabe and others. 
In order to prove the complete equivalence of the 
electrical and mechanical theories, it only 
remains to show that Eqs. (23) and (24) are 
consistent with Eqs. (1) and (2). From Eqs. (23) 
and (24), R’/R=L’/L=(e’/e)*. From the equa- 
tions in §2 for Ci, C2 and C;’ it is evident that 
this ratio is the same as (C,/C,’)?=(1+C,/C2)’, 
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jn agreement with Eq. (1). Agreement with Eq. 
(2) is found from Eqs. (23) and (24) together 
with (32) and (33). It thus becomes evident that, 
so far as thickness vibrations are concerned, 
polarization effects do not invalidate the assump- 
tion that the crystal can be treated as having 
fixed equivalent electrical constants R, L, C and 
C;, with the gap capacitance C, in series. 


(b) Lengthwise vibrations 

§6. Longitudinal vibrations in piezoelectric bars 
have been treated theoretically and experiment- 
ally, by many authors (especially references 1, 
2, 5, 8, 9, 11, 12, 13, 17, 18, 19). The present 
paper extends the theory to a consideration of 
surface charge and gap effects. Since the only 
experiments recorded here have to do with a 
quartz bar vibrating longitudinally in the direc- 
tion of a mechanical axis, the field being parallel 
to an electrical axis (X-cut, Y-waves), the 
terminology applies only to this type, but the 
theory is equally applicable to any sort of piezo- 
electric bar of rectangular section excited by an 
electric field perpendicular to its length. On the 
experimental side most of the sources of error 
mentioned under thickness vibrations may be 
present. 

The theory follows in general the same method 
as for thickness vibrations. (It is assumed that 
the bar is of sufficient breadth for edge effects 
to be neglected.) Using the same notation as 
before, we have*"!-!8 for a bar of dimensions 
e, 1, b parallel to the x, y, z axes, respectively, 
and calling b/=A, 


M=peA/2, W=n%pQbe/2l, 
G’ = r°bec’ /21. 


As in the case of thickness vibrations, the 
logarithmic decrement is A= W/2fM. The driv- 
ing force impressed on the equivalent mass M 
is F= Fy cos wt. The driving force involves the 
effective piezoelectric stress constant e, which in 
the present case is not €1; but 611/511, where sy, is 
an elastic compliance constant, and 61; a piezo- 


™S. Namba and S. Matsumura, Researches of El. tech. 
2. Ministry of Communications, Tokyo, No. 248, April 


Vigoureux, Quartz Resonators and Oscillators (Lon- 


don, 1931 


woes Giebe and E. Blechschmidt, Ann. d. Physik 18, 457 


electric strain constant. The proof of this lies in 
the fact that an electric field parallel to x excites 
the three piezoelectric stresses Y,, X, and Y,, 
the last two of which, through the elastic cross 
constants (Poisson’s ratio) make contributions 
to Y,. From the fundamental piezoelectric and 
elastic relations” as applied to quartz, it follows 
that the total impressed stress causing longi- 
tudinal vibrations parallel to y is (Y,)i=(eu 
—€11(S12/S11) + €14(S14/S11)) which is identically 
equal to 6::£;/s1: (reference 20, p. 817, Eq. (13)). 
This reasoning is correct only when the bar is 
so narrow that vibrations of the modes x, and 
yz have natural frequencies that are high in 
comparison with the y, frequency; otherwise, 
X, and Y, would be met by counter stresses 
that would require modification of the equation 


for (Y,)1. 


We may now write for the driving force, 
letting stand for 61;/S11, 


— F=2ebeE,=2ebe V/e' = (2ebe Vo/e’) cos wt. 


The displacement at the end of the bar is = sin 
(wt— 6), where (neglecting the sign) 


Fo/wZ, = 2ebe Vo/we'Z., 
Z2= W?+(wM—G'/w)?. 


The strain at distance y from the center is 
yy=(ré/l) cos ry/l= Y,/c’, where Yy, is the 
stress at y and c’ is the effective elastic coefficient. 
Since the only stress existing in the bar is Y,, the 
polarization at y (parallel to the x axis) may 
best be expressed as P(y)=61Yy=—d6ic'yy 
= —(m6,,:c'E/l) cos ry/l, where 61; is the piezo- 
electric strain-constant. Since, as will be seen, c’ 
for a quartz bar is at all gap values very nearly 
equal to 1/s11, we may with negligible error 
write 6,:c’ ~6,;/sii=e, as in the case of the 
driving stress. The polarization is distributed 
sinusoidally, being zero at the ends of the bar, 
with a maximum at the center, but it has no 
space variation parallel to the x axis. Hence it 
produces no space charge, but rather polarization 
surface charges o= +P at the boundaries. These 
charges give rise to a field component E, and 
displacement D in the gap, which at any y has 
the value E,,(y) = D(y) =42P(y)e/e’, an equation 


2 W. Voigt, Lehrbuch der Kristallphysik (Leipzig, 1910) 
(or second edition, 1928), p. 901. 
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which may be generally proved for any uniformly 
polarized dielectric plate between parallel elec- 
trodes at zero potential. On substituting the 
value of P(y) we find D(y) = —(41*eet/e’l) cos 
mry/l. D(y) has the same sign as P; like P, it 
has a sinusoidal space variation parallel to the 
length of the bar, but is uniform in gap and 
crystal in the x direction. On the other hand, 
the field E(y) in the dielectric caused by o 
is of opposite sign to P and D. Its value may be 
shown to be 


E(y)=—4P(y)w/e’ = cos ry/l. (35) 


We next find the expression for the current J 
in the resonator, analogous to Eqs. (15) and 
(21), the first term of which is unchanged. For 
the second term we need first the total electric 
flux ¢ over the surface A. This is obviously 


1/2 


g=b Ddy = 8rebet/e’. 


Hence 


1 dy 2ebedé 
Vo cos (wt— 86) 


e’ ran 


V 
=— cos 6), 


where the equivalent electrical impedance is 

Z'={R"+(wL’—1/wC’)*}*. As in Eq. (22), we 

write 7’ =e" /4eb’e?, and have 
R’=Wr', L’=Mr’, 


C’=1/r'G’ (36) 


and for zero gap, 


R=Wr, L=Mr, C=1/rG. (37) 


As in the case of thickness vibrations, so also 
here C,'=kA/4re’. 

It remains to find G for zero gap and G’ for 
gap w. Instead of the two stiffness terms G,; and 
Gz in Eq. (32) we now have only one, designated 
as Gi, since there is but one term in Eq. (35) for 
E. The reacting stress Y,’ due to the depolarizing 
field E, corresponding to X; and Xe, Eq. (26), is 
found from Eq. (35), using the same piezoelectric 
constant ¢€ as in the case of the driving stress. At 
any value of y, 
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4r 
— Y,'(y) =«E(y) = 


Yy=Ciy1. (38) 


The ‘quasi-elastic constant’’ js 
always positive, hence it must be added to the 
pure elastic constant, which for a bar is Young’s 
modulus. For a quartz bar perpendicular to the 
principal (z) axis this constant is 1/s,; (references 
5, 8, 21, 22, 23, 24). The effective elastic constant 
with gap w is therefore 


1 4rew 1 
Sit e’ Si esi? 


At zero gap this reduces to c=cy=1/si1, the pure 
elastic constant. As the gap approaches infinity, 
c’ approaches the value 1/51:+4761:2/ks1:? and 
the frequency increases correspondingly. The 
effective value of the compliance constant, which 
we will call si’, is therefore for infinite gap 
Sir’ in agreement with Giebe and 
Scheibe,> whose deductions, reached by a dif- 
ferent method, are based on the tacit assumption 
of infinite gap. Setting 
—6.9(10-*), k=4.46 (Table I), one finds that 
which is approximately 
one percent of $11. 

The values of the piezo-capacitances C’ and 
C in Eqs. (36) and (37) can now be given in 
detail. By the method described in §4, we find 
G’ = G= 1*be/2lsi1, whence 


C’=8ebel/ne"c’, (40) 


If the vibrating bar, when a gap w is present, 
were strictly equivalent to a bar with zero gap 
in series with a condenser C.=5)/l/4rw, the 
values of R’, L’, C’, R, L and C in Eqs. (36) and 
(37) should satisfy Eqs. (1) and (2). This is at 
once seen to be the case for the ratios R’/R 
and L’/L, but not for C’/C. Eq. (2) is easily 
converted into the form 


C/C’ =(e" /e*)(1+eC/e’C2), (41) 


*1'W. P. Mason, Bell Sys. Tech. J. 13, 405 (1934). 

2 R. B. Wright and D. M. Stuart, Bur. Standards J. 
Research 7, 519 (1931). 

231. Koga, Proc. I. R. E. 24, 510 (1936). 

241, Koga, Proc. I. R. E. 24, 532 (1936). 
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while from Eqs. (36) and (37) we obtain 


e ble’ 8 

The discrepancy lies in the fact that when the 
electric flux in the gap is sinusoidally distributed, 
so that the crystal surface is not equipotential, 
the simple equation for C2, which is assumed in 
Eqs. (41) and (42), no longer holds. The dis- 
crepancy disappears at zero gap, since the 
crystal surfaces are then equipotential. 

These considerations show that in the case of 
a bar driven by a transverse alternating electric 
field it is not permissible to treat the resonator as 
an RLCC, mesh in series with the capacitance of 
the gap, as in Fig. la. It is, however, correctly 
_ represented in terms of the electrical equivalents 
R'L'C’Cy’ given by Eqs. (36); or C’ may be 
derived from Eq. (40) or (42). The dependence 
of frequency upon gap is given by Eqs. (45) and 
(47a). 

§7. Crystal bar metallically coated. If the two 
faces of the bar that are normal to the electric 
field are lightly silvered, thereby causing them 
to remain equipotential even when there is a gap, 
the analysis employed in the last section shows 
that the only equivalent electrical constants to 
be altered are C and C’. This alteration comes 
about through the fact that Y,’ now has a 
constant value throughout the bar, namely the 
value given in Eq. (38) averaged from —1/2 to 
1/2. The corresponding quasi-elastic constant, 
which may be called ce, is derived by the method 
used in deriving Eq. (29) and has the value 
(2=326,,°w/ms17e’. This takes the place of c 
in Eq. (39), from which it-differs by the factor 
8/n*, indicating a slightly lower resonant fre- 


quency, when there is a gap. At infinite gap the 
difference amounts to about one-tenth percent. 
The effective elastic constant is c’=co+ce, and 
this is the value to use in computing G’ and C’ 
for a silvered plate with gap. Cs assumes at all 
gaps the simple form C.,=6//4rw. In this case, 
and only then, can the resonator be regarded as 
equivalent to Fig. 1a, C2 representing the geo- 
metrical capacitance of the gap. 

In all the above expressions for capacitance it 
must, of course, be remembered that the formula 
for infinite parallel plates ceases to be even ap- 
proximately correct when the gap length w 
approaches the crystal width 6 in value. Reliable 
experimental tests must, therefore, be confined 
to relatively short gaps. 


II. Various Cuts AND THEIR FUNDAMENTAL 
CONSTANTS 


List of cuts employed 


§8. In Table I the first column indicates the 
cuts for which experimental data were obtained. 
“x” and “‘y” signify plates perpendicular to the 
x and y oe respectively. The AC- and BT-cuts 
are those described by Lack, Willard and Fair ;?5 
similar investigations were published almost 
simultaneously by other writers.**-* For a right 
quartz crystal the last two cuts may be con- 
sidered as obtained by rotating a Y-cut plate 

%F. R. Lack, G. W. Willard and I. E. Fair, Bell Sys. 
Tech. J. 13, 453 (1934). 

6 R. Bechmann, Naturwiss. 21, 752 (1933); Hochfreq. 
tech. u. Elek. akus. 44, 145 (1934). 

27 1, Koga and N. Takagi, J. Inst. Elec. Eng. of Japan 53, 
940, 1141 (1933). 


28 J, Koga, Rep. of Radio Res. in Japan 4, 61 (1934); 
— No. 33, Comm. I, U. R. S. 1. Gen. Assembly, London 
934. 


ot Straubel, Physik. Zeits. 35, 179 (1934); see also 
Zeits. f. tech. Physik 15, 607 (1934). 


TABLE I. Data for various cuts. For the values of co, see Table III. The magnitudes of the : on eages constants ¢; 


and cz, Eq. (34), may be derived from the last two columns. For the quasi-elastic constants 


the bar, see §§6 and 7. 


TYPE OF VIBRATION Type oF STRAIN’ k 


4ré 32é 
€ k wk 


thickness compressional Xz 4.42 
thickness shear 4.42 


lengthwise compressional 4.46 


«104 108 108 
é1.=—5.1 Cu 74 60 


é1.=—5.1 cw 60 


= —3.16 Coe’ 23 
€26 = —2.86 Coe’ 5 
611/s11 = —5.33 1/siy ome 
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about the x axis through an angle 6, this angle 
being positive when the rotation is clockwise on 
sighting along the positive direction (Voigt’s 
definition, reference 20, pp. 750, 859) of the 
x axis. For a left quartz the sign of @ is to be 
reversed. The value of this angle found by Lack, 
Willard and Fair for a right quartz is for the 
AC-cut 31°, and for the BT-cut, —49°. 

In the third column of Table I is indicated 
the type of strain (and therefore of stress) for 
each cut. Primed symbols indicate rotated axes. 


Elastic constants 


§9. For thickness vibrations, the proper elastic 
constant for each cut, which in generalized form 
was denoted by co in §4, is that one of the c’s (in 
Voigt’s notation) which expresses the stress- 
strain relation characteristic of the type of 
vibration.”: *! In quartz, the difference between 
the adiabatic and isothermal values is less than 
that between isothermal values as reported by 
different observers. The reason why c should be 
employed rather than the reciprocal of a com- 
pliance constant s is, that while there is but one 
strain, there is a system of lateral stresses (con- 
straints) to which the absence of other strains 
is due. This statement is of course rigorously 
true only of infinite plates. 

For lengthwise vibrations of a narrow bar 
parallel to the y axis, the expression for Young’s 
modulus is 1/511, as was stated in §6. 

The specific form assumed by the generalized 
constant Co for each of the cuts discussed in §19 
is indicated in the fourth column of Table I. 
Numerical values are given in §19. For the 
derivation of ces’ see reference 25, p. 463. The 
primes indicate oblique cuts with transformed 
axes. 

In carrying out the final computations it was 
deemed wisest to derive the effective values c of 
the elastic constants from the measured fre- 
quencies, using the formula c=4pf’e? (or 4pf?/* 
for the bar) where the density p=2.65. These 
values are shown in Table III for each crystal. 
From them, with the aid of Eqs. (34) and (39), 
the pure elastic constants cy have been calculated 
(see Table III). For comparison it may be noted 
that Voigt’s values (adiabatic) for ¢11, Ces and 
1/s\, are, respectively, 8.54, 3.91 and 7.72, all 
times 10" dynes/cm?. 
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The dielectric constants 


§10. Whenever the dielectric constant of a 
piezoelectric crystal is measured by the usual 
method of placing a plate of the material in an 
electric field, the observed value is always too 
great (reference 20, p. 917; reference 1, second 
citation, p. 749), owing to the deformation and 
consequent polarization caused by the electric 
field. The pure dielectric constant could be 
directly measured only if during the measure- 
ment the plate were kept unstrained, which is 
practically impossible. 

On applying to Voigt’s Eq. (287) (reference 20, 
p- 917) the subscripts for an X-cut in quartz, 
there results for the polarization produced by a 
field E, Ei, for 
static conditions with the plate devoid of external 
stresses. m1. is the dielectric susceptibility . 
parallel to the x axis. With thickness vibrations 
in an infinite plate each unit of volume is 
subject to external stresses due to the lateral 
constraints which neutralize the last two terms. 
The strain corresponding to the second term, 
which is of the type x,, is replaced by the vibra- 
tional strain from which the piezo-current i, 
Eq. (15), is found, so that the only pure dielectric 
polarization caused by the impressed field is 
from which it follows that k=1+477, 
the value for an unstrained crystal. The best 
value of the observed dielectric constant of 
quartz perpepdicular to the optic axis (crystal 
free from constraints) is 4.5.° From Voigt’s 
equation, and setting 1611+ 2612+ 414614 = 0.0068, 
we have whence k=1 
+4nn=4.42 for the unstrained crystal. Since the 
susceptibility » is the same in all directions per- 
pendicular to the optic axis, it follows that 
k=4.42 is the proper value for both X- and Y-cut 
plates (Table I). In the case of the AC- and BT- 
cuts it is necessary to compute k from the 
observed values parallel and perpendicular to 
the optic axis, namely 4.6 and 4.5, respectively.” 
The results are recorded in Table I. The values 
of k in different directions are so nearly alike 
that refraction of electric lines at the surface of 
the quartz (when there is a gap) may be ignored. 

In the case of lengthwise vibrations, the 
stresses corresponding to the terms €1:61: and 


* R. B. Sosman, The Properties of Silica (New York, 
1927), p. 515. 
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€:4614 are not neutralized by external constraints, 
whence it is found that the effective dielectric 
constant is 4.46. 

Although these distinctions are too small to 
be easily noticeable in quartz, it seems worth 
while to call attention to them, since they may 
easily play a more important part in resonators 
employing other crystals. 


Piezoelectric constants 


§11. The values adopted for ¢:; and 6,; are 
taken from the International Critical Tables, 
Vol. VI: €11.= —5.1(104), —6.9(10-). The 
sign applies to left quartz; for right quartz it 
should be reversed. The distinction is, of course, 
of no importance in the present vibrational 
equations. For s,; in Table I, Voigt’s value, with 
the adiabatic correction, is used : = 12.95(10-"*) 
cm*/dyne. The values of és’ for the AC- and 
BT-cuts were computed by Voigt’s method 
(reference 20, p. 840), modifying his equations 
for rotation about the x axis.*! 


III. GRAPHICAL TREATMENT OF THE 
RESONATOR 


Resonance circle diagram 


§12. The relations between various parameters, 
as well as the contrast between the electrical and 
the mechanical theories, can best be understood 
by means of the resonance circle diagram, Fig. 4. 
The treatment previously described® will be 
employed, according to which the same diagram 
serves for both admittances and impedances. We 
shall first explain the construction of the diagram 
representing Fig. 1a, according to the electrical 
theory, and then interpret it in terms of the me- 
chanical theory. When there is no gap, the circle 
is based on the electrical constants R, L and C, 
Fig. 1a, which were shown in Eqs. (24) to be pro- 
portional to the mechanical constants W, M and 
1/G. It therefore represents in electrical terms the 
state of vibration of the resonator. In order to 
adapt the diagram to both impedances and admit- 
tances, the two unit vectors (scale-values) s, and 
5, are introduced, representing, respectively, the 


“In the 1910 edition of Voigt the wrong sign occurs 

fore the second parenthesis in the equation for e2.', 
p. 841. This has been corrected in the 1928 edition. 
®W.G. Cady, Proc. Am. Ac. 68, 383 (1933). 
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Fic. 4. Resonance circle er Origin is at F’ when 
gap is present, approaching F as the gap diminishes to 
zero. 


number of reciprocal ohms (‘‘mhos’’), and the 
number of ohms, per unit distance on the 
diagram. If pp=AC is the radius of the circle in 


‘terms of unit distance, then s,-2p9.=1/R, and 


the length of the line from A to any point on the 
circle, multiplied by s,, represents to scale the 
admittance Y of the RLC branch at the fre- 
quency corresponding to that point. For example, 
at point 7, Fig. 4, the RLC admittance is 
Y=s,-A7, having components g=s,:AD and 
b=s,-D7. 

From Eqs. (9), (13) and (24), it is evident 
that any line from A to a point on the circular 
locus, as for example A7, with proper choice of 
scale values (Table II), also represents the 
mechanical admittance Y,=s,-A7, as well as 
the maximum velocity of motion, vo=s,-A7. 
S$, is expressed in cm/sec. per unit of V per scale 
division. Velocity, mechanical admittance and 
equivalent electrical admittance, all have the 
same circular locus and, for any frequency, the 
same phase angle 0, defined by Eq. (12). Positive 
values of 6 lie above the horizontal, indicating 
an inductive equivalent circuit for RLC. 

From Eq. (13) it is seen that the amplitude of 
vibration &) has also very approximately a 
circular locus. The maximum amplitude does not 
come at the frequency corresponding to point B, 
however, but at a very slightly higher frequency, 
the angular velocity being (G/M—W?/2M?)} 
= (wo?—2a?)!, where wo is the angular velocity 
at B and a is the damping constant (reference 12, 
p. 5). The scale value for & is that for v» divided 
by w. 
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The frequency distribution around the circle 
depends on R, L and C (not at all on C;), and 
runs counterclockwise for admittances, as indi- 
cated by the arrow in Fig. 4, and clockwise for 
impedances. In the case of the admittance 
diagram for RLC the frequency is zero at A, 
increasing over the lower half of the circle to the 
resonance value fo at B, and finally reaching the 
value infinity at A. 

The vertical line through B serves as a con- 
venient frequency scale, the scale being very 
approximately linear over the resonance region 
(reference 33, p. 135; reference 32, p. 390). It is 
only necessary to prolong the line from A through 
the point on the circle for which the frequency 
f is desired until it meets the vertical line through 
B at some point S. Then to a close approximation 
n=a-BS, where n=f—fy and o=R/(8rpoL). 
BS is positive for points above the AB axis. As 
an illustration, the construction has been carried 
out in Fig. 4 for point 3. For points on the left 
half of the circle the distance BS becomes incon- 
veniently large. In this case an alternative for- 
mula may be used (reference 32, p. 391): 
n=a'/AV, where o’=poR/2xL and V is the 
point on the vertical axis through A where the 
prolongation of the line from B through the 
point on the circle meets this axis. This case is 
also illustrated in Fig. 4 for point 9, V being the 
corresponding point on the frequency scale, and 
n=a'/AV being the frequency-difference fz —f». 
The AV scale is not linear, and like the BS scale 
it is accurate only over the resonance region. 
Because of the very low decrement of the 
resonator, all of the circle except a small region 
close to A corresponds to frequencies differing 
but little from fo. fo is that frequency for which 
the reactance of the RLC branch vanishes 
(4°f,?=1/LC); it is therefore the frequency of 
maximum velocity of motion under constant 
driving force (Eq. 13). 


The RLCC, mesh 


§13. In order to represent the admittance Yj’, 
of the RLCC, mesh at the frequency correspond- 
ing to any point P on the circle, we lay off 
AF=-—w(C,/s, along the vertical axis. The 
upward direction is chosen in order that, when 
F is taken as the origin of vectors, capacitive 


3% E. Mallett, Telegraphy and Telephony (London, 1929), 


WALTER G. CADY 


values will come below the horizontal. The fre- 
quency varies so little over the resonance range 
that the points F and F’ may be regarded as 
fixed. The vector sum FP of FA and AP repre- 
sents Y,’, P being any point on the circle. For 
example, at frequency fo, P coincides with B 
and s,-F B= Y,’; this shows that at the frequency 
for maximum velocity the electrical admittance 
of the crystal is slightly capacitive. At all fre- 
quencies above the horizontal line through F 
the crystal admittance is inductive. At points 1 
and 2 the crystal behaves as a pure resistance. 
The admittance at 1 is large, and the frequency 
is very slightly above fo. At 2, F2 is small and 
the condition is that of parallel resonance; the 
admittance F2 is for most resonators much 
smaller, even, than is shown in Fig. 4. This is 
because the distance AF has been greatly exag- 
gerated for clearness: for a good resonator the 
ratio AF :AB=(CiR is of the order 1 : 100. 
The frequencies for maximum and minimum 
admittance Y,,’ and Y,,’ fall at the points 3 and 4 
(vector from F through C). At 4 is the condition 
of ‘“‘antiresonance.”’ The analytic expressions for 
the various vectors are cumbersome and _ need 
not be given. Dye® gives an extensive analytic 
treatment, while the general equation for the 
admittance of a resonator with gap has been 
derived by Usui.* 


Crystal with gap 


$14. The method of deriving an admittance 
diagram for R’L’C’C;’, Fig. 1b, analogous to 
that for RLCC, consists first in inverting the 
latter, which may for brevity be called the s, 
diagram, into an impedance, or s., diagram, 
keeping F as origin. The origin is next shifted 
to F’, the location of which is determined by the 
equation s,- AF’ = —1/wC;,’ or s,-FF’= —1/wC, 
where s., the impedance scale value, is 1/s,(AF)? 
(reference 32, pp. 389, 396). We thus obtain the 
impedance diagram for RLCC,C2. One more 
inversion, this time with respect to F’, is now 
performed, giving the desired R’L’C’C,’ ad- 
mittance diagram (the s,’ diagram). F’ becomes 
the origin of admittance vectors, and s,'/AF’= 
—w(C;’. The value of R’ is expressed by an equa- 
tion analogous to that for the s, diagram: 


4 R. Usui, J. Inst. Elec. Eng. of Japan 54, 201 (1934) 
(in Japanese; summary in English, ibid., p. 21). 
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TABLE II. Resonance diagram notation. 


NETWORK ORIGIN VECTOR 


ELECTRICAL 


SCALE VALUES 
MECHAN. ADMITTANCE VELocITY 


RLC A Y 


RLCC; F 


RLCCi F 
RLCC\C: F’ 


RLCC,C2 


RLCC;C2C3 F;’ Y;’ 


RLCC\C2C5 


Sy =1/2poR 


e 


=1/(AF3’)?s, 


_ e 
SepoA?R 


s:=1/(AF)*s, 


s,/-2p0=1/R’. The quantities involved are sum- 
marized in Table II. . 

Consider, for example, point 3, Fig. 4, at a 
frequency f; slightly less than fo. We have 
Y=s,-A3, Y;'=s,:F3. According to the prin- 
ciple of inversion, Z,’, the reciprocal of Yj’, is 
obtained from 4, the point inverse to 3 with 
respect to F, that is, Z;/=s,-F4. The vector 
sum of F’F and F4 is F’4, or Z.’=s,-F’4. By 
inverting 4 with respect to F’ we obtain point 9, 
from which the admittance, at the same fre- 
quency, of the entire network RLCC,C2, may be 
found: 

The remarks that were made concerning points 
1, 2, 3, 4 and 9 in $13 hold, when there is a gap, 
with respect to points 5, 6, 7, 8 and 10. If the gap 
is so great that AF’>AC, the reactance of the 
resonator is capacitive at all frequencies. 

The circle diagram with F’ as origin thus 
represents the performance of the R’L’C’C,’ 
network, Fig. 1b, and from Eqs. (23) it is evident 
that it represents also the state of mechanical 
vibration, just as was the case with zero gap. 
In particular, the frequency at B, which may be 
called fo’, is that at which the velocity v’, and 
also the piezo-current 7’, is a maximum. Although 
it is not recognizable as a critical frequency by 
direct electrical observation, still fo’ is the 
resonant frequency that is most accurately 
characteristic of the piezoelectric behavior of the 
resonator, and as such is of fundamental im- 
portance. Being practically identical with the 
frequency of maximum amplitude £, it is 
the frequency observed in experiments using the 


interferometer or microscope. The greater the 
gap, the more is the value of fo’ in excess of fo, 
the value for zero gap. This difference, according 
to the mechanical theory, is due to the depend- 
ence of the stiffness G’, Eqs. (14) and (32), upon 
the gap. It is this recognition of the variation of 
the effective mechanical stiffness with gap, 
rather than the explicit use of the gap capaci- 
tance, that distinguishes the mechanical from 
the electrical theory. The advantage of the 
mechanical theory lies in the better description 
that it gives of the physical phenomena, as well 
as in the greater ease with which the vibrational 
state may be determined at any frequency, and 
therefrom the potential difference across the 
crystal and the danger of breakage (§15). 

As the gap diminishes, F’ approaches F, and 
the resonator with zero gap becomes a special 
limiting case. Since the R’/L’ ratio, and hence 
the decrement, are independent of the gap, it 
follows that the same frequency scale value o 
($13) holds at all gaps.* The quantities de- 
pendent on the gap are the position of the 
origin F’, the admittance scale value s,, the fre- 
quency fo at B, and the distribution of fre- 
quencies around the circle. With a given gap, 
if a point on the circle is experimentally de- 
termined at a frequency f, the frequency fy’ 
may be found immediately by graphical con- 
struction, as was illustrated for points 3 and 9, 
Fig. 4, in §13. By performing two obvious in- 


% Watanabe (reference 3, second citation, p. 710), finds 
experimentally an increase in decrement at very small gaps, 
due apparently to air friction. 
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versions with respect to F’ and F, it is easy to 
locate on the s, circle the frequency fo’ at point 
B when there is a gap, or indeed the frequency at 
any point whatever on the s,’ circle. If the 
frequency f,’ for maximum admittance with 
gap has been observed (point 7), the graphical 
method shows at a stroke (by drawing F8) that 
the same frequency comes at point 10 on the sy 
circle, which, as it depends on the fundamental 
constants, may be regarded as the fiducial 
circle. In general, the greater the gap, the more 
are the frequency values shifted clockwise around 
the circle. 


Phase diagram 

§15. The relative phases and magnitudes of 
the various periodic quantities involved in the 
resonator are shown in Fig. 5. As may easily be 
verified, they are derived from Fig. 4, for the 
frequency at point 7 on the admittance diagram 
with gap, origin at F’. For example, since F’7 
represents the total admittance Y,’, it is evident 
that the total current J (Fig. 1) leads V by the 
angle 6F’8. The various vectors may be regarded 
as either maximum or r.m.s. values. 

It is necessary first to explain the voltages V. 
and V;. V; is the component of potential 
difference between the crystal surfaces (Fig. 3) 
due to the state of strain, on the assumption of 
zero p.d. between electrodes. It is therefore equal 
and opposite to the corresponding drop across 


_Fic. 5. Phase diagram for resonator with gap, based on 
rig. 4. hag save space, the vectors i and 7, are drawn to 
scale. 
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the gap, which by Eq. (20) is 8xetw/e’. V. is 
the actual p.d. between the crystal surfaces and 
is, therefore, the vector sum of V; and Vj, as 
shown in Fig. 5, where Vi=eV/e, is the primary 
voltage impressed on the crystal (Fig. 2). A 
general equation for V. in terms of equivalent 
electrical constants is easily derived from the 
expression above for V;, with the aid of Eqs. (4), 
(9) and (11): 


4rwY’ e 
Ve= sin (wt—6)+— cos ut}. 


w 


The phase angles, derived from Fig. 4, are as 
follows, the notation ‘i,’ : V=90°” signifying 
that 7,’ leads V by 90°, etc.: I: V= Z6F’8; 
I: Ve=90°; 4,’ : V=90°; V= ZBA7; V.: 1 
= Z8F1; V.:i and V.: v= ZBA10; v: €=90°; 
&£=180°; 1, V.=90°. 

Fig. 5 shows the current components for both 
Fig. 1a and b. Thus J is the resultant of 7 and i, 
and also of 7’ and i’. 

According to the mechanical theory, V, is the 
driving voltage at point 7, Fig. 4, and &, V; and 
V. are large because of approximate mechanical 
resonance, v being nearly in phase with V. 
According to the electrical theory, the crystal-gap 
combination is nearly in series resonance, the 
crystal reactance being strongly inductive, as is 
shown by the angle 10AB, Fig. 4. Point 10 is 
found by inverting 7 into 8 with respect to F’, 
then inverting 8 into 10 with respect to F; it is 
the point on the RLC admittance circle at which 
the frequency is the same as at 7 on the R’L’C’C,’ 
admittance circle. Thus according to the elec- 


trical theory the crystal itself is far from reso- 


nance, but the amplitude is large because V,, 
which is now regarded as the driving voltage, is 
large. This is in striking contrast with the 
mechanical theory, which says that the drop 
across the crystal is great because the amplitude 
of vibration is large; we have an apparent inter- 
change of cause and effect! The reason lies in 
the fact that the mechanical theory regards the 
effective stiffness as dependent on the gap, as 
has already been explained. The relations would 
appear still more striking if it were feasible to 
draw vector diagrams for a normal crystal, in 
which the ratio wC, : 1/R and hence AF : AB is 
of the order 1 : 100 or less. 
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The production of very large voltages across 
the gap when the resonator is in the state of 
optimum gap resonance has been discussed by 
Dye? and the author.* 


IV. EQUATIONS FOR CHANGE IN FREQUENCY 
WITH GAP 


Derivation of coefficient U 


§16. Since the frequency of greatest theoretical 
importance is fo’ (§14), an equation will first be 
derived for Af=fo’—fo, where fo (§12) is the 
limiting value of fo’ at zero gap. From the 
definitions of fo and fo’, and the relation 
4e*f, o =c'/p, 


For thickness vibrations, from Eqs. (33), (34) 
and (43), it is easily found that 


Af/fo= 
or, remembering that e’=e+kw, 
Af/fo=16¢?w/mce’ (44) 


where c is the value of c’ at zero gap. In the 
limit, when w= this becomes 16¢2/zck. 

For lengthwise vibrations, when the surfaces 
of the bar are not made equipotential by a 
metallic deposit, it is found by the above 
procedure, using Eqs. (36), (37), (39) and (40), 
that 


Af 2nréw w 
— = — = 276, (45) 
fo co 


An equation of the same form as (44) is found 
in the case of lengthwise vibrations when the 
surfaces of the bar are made conductive. The 


_ same equations are employed as in the derivation 


of Eq. (45), only in this case the expression for ce 
given in §7 takes the place of that for c; in Eq. 
(39). There results the equation 


Af 16e w 166,:° 
—. (46) 


This expression is smaller than that in Eq. (45) 
by the factor 8/7. 

The last three equations may be written in 
the general form Af/fyo= Uw/e’, the coefficient U 
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in each case being a function of piezoelectric 
and elastic constants alone. It was the object of 
the experiments to compare this theoretical with 
the observed value of U. 

When the crystal surfaces are equipotential, 
the effect of gap on frequency can also be ex- 
pressed in terms of the gap capacitance C, 
(reference 2, p. 424; reference 18, p. 51). Thus 
from Eq. (44) we find for thickness vibrations, 
with the aid of Eqs. (24) and (33), together with 
the equation w/e’ =C,/k(Ci+C2), 


Af/fo=3C/(Ci+C2). (47) 


For infinite gap this becomes C/2C,. 

The same equation can similarly be derived 
for lengthwise vibrations (crystal silvered) from 
Eq. (46). In the case of lengthwise vibrations, 
when the bar is not metallically coated, it is 
found from Eq. (45) that 


Af/fo= (x? (47a) 


The piezo-capacitance C can be measured 
electrically without reference to ¢ or 38 
When this is done, Eq. (47) can be made to 
serve as a check on Eqs. (44), (45) and (46), and 
indeed it becomes a means of determining the 
piezoelectric coefficient « at high frequency. No 
such observations were made with the crystals 
discussed in the present paper. 

The theoretical relative change in frequency 
Af/fo, as the gap increases from zero to infinity, 
is easily found from the equations above by 
substituting for w/e’ its limiting value 1/k. 
It is approximately for the Y-cut 7.2; B7-cut 
1.3; AC-cut 3.9; X-cut, thickness vibrations, 
3.39; X-cut bar, unsilvered 5.1, silvered 4.2, all 
multiplied by 10-*. The values of 1/k for the 
above cuts (see Table I) are, respectively, 0.226, 
0.221, 0.224, 0.226, 0.224 and 0.224. 


V. OBSERVATIONS AND RESULTs*® 
Apparatus 


§17. With all thickness-vibration crystals a 
variable-frequency oscillator, the frequency of 
which could be controlled and read to one part 


36 The experimental work was carried out by Mr. John F. 
Miiller. The writer wishes to express his thanks to the 
Bell Telephone Laboratories for the loan of certain crystals 
and for other courtesies. 


| 
| | 
) 


254 


in a million, was connected to the points a, b in 
Fig. 6, so that the resonator was excited by the 
voltage drop across the resistance Roy of 600 
ohms. Harmonic frequencies were removed by 
filters. For each observation the voltage was 
kept constant while the frequency fn,’ was 
recorded at which the galvanometer G connected 
to a thermoelement indicated maximum current, 
(usually less than 0.5 ma), and hence maximum 
admittance of the RLCC,C2C;R; network. The 
heater resistance R, had a value of 500 ohms. 
C; represents the stray capacitance of lead wires 
(Table III). Each crystal was carefully ground 
to uniform thickness and lay upon the lower 
electrode in a special mounting. The gap was 
regulated by a micrometer screw attached to 
the upper electrode. The faces of crystal X 1 
were plane-parallel to within 4 microns, the 
other thickness-vibration crystals to within 1 
micron. Tests with polarized light revealed no 
traces of twinning. . 

For crystal Q 30 (lengthwise vibrations) a dif- 
ferent oscillator was used, and the difference 
Afm’ between the frequency at each gap and 
that at the smallest gap was observed. The same 
measuring circuit, Fig. 6, was used. This crystal 
rested horizontally upon the flat end of a round 
brass rod 1.5 mm in diameter which extended 
upward through the center of the lower electrode 
for a very short distance. It was clamped in 
position by a second brass rod 0.8 mm in diameter 
extending downward through the center of the 
upper electrode. The gap was varied by means 
of a micrometer screw which controlled the upper 
electrode. Only the fundamental lengthwise 
compressional mode was observed. Previous tests 
by dust-patterns had shown this to be the only 


w 


element. 
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perceptible mode. The resonant frequency, 
66,500 cycles per second, was not appreciably 
altered by silvering. The designations Q 30 and 
Q 30s in Table III refer to this crystal in the 
unsilvered and silvered conditions, respectively, 
In all cases the electrode faces were identical 
in size with those of the crystal. All the upper 
electrodes were made accurately plane, as were 
also the lower electrodes of crystals X 1 and Q 30, 
The lower electrodes of the remaining crystals 
were slightly concave, to a depth of not over 
0.025 mm, in order that the central portions 
might be free to vibrate. Every precaution was 
taken to minimize stray capacitance. Approx- 
imate values of the equivalent resistances R had 
previously been determined (Table III). Values 
of L were computed from Eqs. (14) and (24). 


Derivation of Af from observed frequencies 


$18. Because of the presence of C; and R, the 
observed frequency f,’ is not quite the same as 
fo’ in Eq. (43). If the resistance R, is negligible, 
fo’ is best derived from f,,’ graphically, except . 
with small gaps, when approximate formulas 
may be used. The graphical method is indicated 
in Fig. 7, in which, as in Fig. 4, F is the origin 
of vectors for the RLCC,; mesh (zero gap), F; 
that for RLCC,C; (zero gap, but including the 
stray capacitance C;) and F’ the origin for the 
RLCC,C:2 network, or its equivalent R’L'C’C,’, 


Fic. 6. Measuring circuit with resonator and thermo- Fic. 7. Circle diagram for reduction of frequency obser- 
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Fig. 1. Since C; is in parallel with C,’, the origin 
for R'L’C’C1'C; is at F;’, where s,’-F’F3’= —wC; 
and 

AF3’= — C;). (48) 


Similarly, the point F; is determined by the 
equation 
sy: AF;= —w(Ci+Cs3). (49) 


At the frequency corresponding to any point P 
on the s,’ circle, the admittance and impedance 
for RLCC,C2C; are, respectively, 
and Z;'=s,'-F3'Ps2, P2 being the point inverse to 
P with respect to F;’. 

If there is no series resistance Rx to consider, 
the frequency f,,’ for maximum admittance falls 
at point P;, Fig. 7, whence fo’=fn,’+o-BS; 
(§14), fo’ being the frequency at B. For zero gap 
(origin at F;), the corresponding equation is 
fo=fm+o-BSo where So is the limiting position 
assumed by S; when the gap is zero. 

The following expression, based on Fig. 7, 
makes it possible, when R, is negligible, to find 
BS;/po and BSo/po without actually performing 
the graphical construction: BS;/po=tan (BASs) 
=2 tan {3 (AF;'/po)}. For each gap, in- 
cluding zero gap, AF;’ is to be found from Eq. 
(48). 

When, as in the present case, the series 
resistance R, is large enough to have an appre- 
ciable effect on the resonance frequency, the 
graphical method offers a great advantage over 
the laborious algebraic process. The first step 
consists in inverting to the impedance diagram 
for R’L'C’C,’C3, the scale value being s,.’, from 
Table II. This quantity, like s,’, has a different 
value for each gap. For any gap, when F;’, 
Fig. 7, has been located, the distance O’F;’ is 
laid off horizontally to the left, by making 
s,-O’F;'=R,, or (see Eq. (48) and Table IT) 


O’F,’ wR,(Cy’ C3)AF;’. (50) 


O’ is then the origin of vectors for the entire 
circuit, including Ry. The locus of O’, as the gap 
varies, can be shown to be a parabola passing 
through A. The curve is most easily drawn by 
calculating O’F;’ for two or three assumed 
values of w. At zero gap, O’F;’ becomes OF3. 
When, for any gap, F;’ and O’ have been 
located, the line O’P; drawn through C repre- 
sents the maximum admittance Y,,’ of the entire 
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circuit, and its inverse, O’Ps, represents the 
minimum impedance at the observed frequency 
fm’, since s,’-O’P,=Z,,. The impedance, at this 
frequency, of the RLCC,C2C; network is s,’- 
and the admittance is s,’-F3;’P. From Fig. 7 it 
is seen that the frequency at P (which includes 
the effect of Rs) is lower than that at P;, the 
point for maximum admittance when R, is 
negligibly small. The value of fo’ ($16) when R, 
is present is 

fo’ =fm' +o0-BS’, (51) 


where as before c=R/8rpoL. At zero gap, the 
same procedure is followed for finding fo=fn 
+o-BS, using O and F; in place of O’ and F;’. 
Since AF;, and also, for small gaps, AF;’, is 
usually small in comparison with the radius po, 
it is best at this end of the range to use the 
geometrical relations indicated in the diagram 
merely as a means for calculating BS’ for each 
gap, including the zero value. The equation is 


BS’/2po=tan gs=tan (¢g1— 


=tan (¢gi—¢2/2). (52) 


The angles ¢; and ¢2 are easily expressed in terms 
of known quantities. This method was pursued 
in deriving the frequency differences Af, with the 
exception of those for larger gaps, when the 
graphical method became necessary. 

Finally, from the equations given above, the 
value of Af for each gap (§16) is expressed, from 
Eq. (51), as 


Af=fo' (53) 


where f,, and BS are the limiting values at zero 
gap. In practice, since it was not feasible to make 
observations strictly at zero gap, the values of 
fo’ were plotted against w/e’, whence fp was found 
by extrapolation. 

The magnitude of the o- BS’ correction ranged 
from 3~/sec. for crystal Q 30 to 151~/sec. for 
X 1. If this correction had not been applied, the 
value of Af would have been lower by about 16 
percent for X 2, 6 percent for X 1, and for the 
remaining crystals from 1.3 to 0.15 percent. 


Results of observations 


§19. The results are shown in condensed form 
in Table III and Fig. 8. The ‘‘calculated”’ 
(theoretical) values of U are the coefficients of 
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TaBLE III. Crystals tested and experimental results. The first three crystals are Y-cut, the next two are cut according 
to Lack, Willard and Fair (§8). X 1 and X 2 are X-cut plates, Q 30 an X-cut bar for lengthwise vibrations, Q 30s the 
same silvered on both sides. For the BT and AC-cuts the x, y, z dimensions should be primed. Crystals Y 3 and X 1 
are circular. Line 4, effective elastic constant at zero gap, according to §9. Line 9, correction f,,’—f’, §18, in cycles per 
sec., for maximum gap used. Line 11, wave constant in meters elmg. wave-length per mm in direction of vibration. 
Line 14, observed U for very small gap and for maximum gap. Line 16, percent excess of calculated over observed U. 

The average values of the pure elastic constants co (see §9) are: Y-cut, 4.15; X-cut (thickness vibrations), 8.81; BT-cut, 
6.92; AC-cut, 2.88; X-cut bar (lengthwise vibrations) 7.66, all times 10" dyne/cm?. For the bar, s;; =1/7.66 (10") =1.30 
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w/e’ in Eqs. (44), (45) and (46), taking « from 
Table I and c¢ from Table III. ‘‘Observed” 
values of U are the ratios of observed Af/fo to 
w/e’. According to theory, U should have a 
constant value for each cut, represented by a 
straight line of constant slope through the 
origin in Fig. 8. This expectation is approximately 
realized with crystals Y 1, X 2 and Q 30, while 
the characteristic lines for the remainder show 
a more or less decided curvature, indicating a 
decrease in U with increasing gap. Nearly always 
the observed values of U are too small, especially 
with the X-cut thickness-vibration crystals X 1 
and X 2. In all cases, with the exception of the 
bar Q30, the original curves showed, over 
certain gap ranges, effects of coupling within the 
crystal, or of air resonance, or both. The air 
resonance effect was very pronounced with the 
X-cut plates, as would be expected. In Fig. 8 the 
detailed curve is shown for crystal Y 1, illus- 
trating strong mechanical coupling within the 
crystal. 

In constructing Fig. 8, smooth lines were 
drawn through the plotted points, to correct as 
far as possible for these effects. There is no 
obvious reason why the lines so drawn should 
not be characteristic of the primary mode. On 


the original diagrams these lines were extended 
until they cut the Af/f axis. In each case, from 
the position of the intercept the basic frequency, 
fo, was determined, from which the corrected 
value of Af was found. It is most unfortunate 
that observations could not be made with the 
gap actually equal to zero, since the process of 
extrapolation led to some uncertainty in the 
precise determination of fo, and hence of Af. 

By continuing the characteristic lines in Fig. 8 
until they meet the ordinates for infinite gap 
(abscissa value 1/k), one can gain an idea of how 
far the greatest possible increase in frequency 
with gap for the various crystals would fall 
short of the theoretical values for the various 
cuts given at the end of §16. This has been done 
in Fig. 8 for crystal Y 1 (Point P). 

The relatively small values of U, even at small 
gaps, may have been due to one or more of the 
following causes: too large an assumed value for 
the piezoelectric constant in the theoretical 
value of U; traces of twinning; and the fact that 
the plates were not of infinite lateral extent. 
This last consideration must have tended, at any 
given gap, to make the frequency less than that 
for an infinite plate, in three ways, as follows: 
First by allowing some lines of electric intensity 


BT AC x1 Q 30 Q 30s 
3.8 3.8 0.286 0.152 0.152 
af 0.335 0.335 4.5 4.04 4.04 
tT 3.8 3.8 4.5 1.40 1.40 
T 6.917 2.890 8.670 7.662 7.662 
200 200 700 2470 2470 
Tt 3.01 0.777 0.170 46.8 46.8 
: 17.2 17.0 21.7 1 13.2 13.2 
8.0 8.0 8.0 6.5 6.5 
M4 4 11 151 91 3 3 
i 7 761,974 492,550 1,000,422 727,598 66,500 66,500 
7 4 118 105 103 112 112 
hi 1.26 1.45 0.85 0.98 0.084 0.25 0.48 0.47 
l. ; 5.27 3.56 2.93 0.76 0.78 0.53 4.15 3.35 
0.0268 0.0238 0.0278 0.0060 0.00875 0.00318 0.0200 0.0158 
0.0268 0.0182 0.0139 0.0038 0.00615 0.00318 0.0200° 0.0158 
a 0.0324 0.0305 0.0327 0.00594 0.0153 0.0147 0.0229 0.0186 
17 15 1 3 43 78 13 15 
» 17 40 58 +36 29 60 78 13 15 
rh 
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to turn back around the edges of the plate, thus 
reducing the depolarization ; second by increasing 
the effective capacitance of the gap, which has 
the same electrical effect as a reduction in gap 
thickness. These two effects are negligible at 
zero gap, but they grow in importance as the gap 
increases, thus offering a partial, if not complete, 
explanation of the curvature of the characteristic 
lines. The fact that the observed curvatures bear 
no consistent relation to the dimensions of the 
crystal plates indicates that this explanation is 
not of itself sufficient. The third consequence of 
the finite dimensions of the crystal is that, as 
the edges of the plate are approached, the lateral 
constraints mentioned in §9 diminish, resulting 
in a diminution in stiffness and hence in the wave 
velocity. The effect is the same as if the thickness 
of the plate increased slightly to a maximum at 
the edges. The natural frequency is thereby 
lowered, and owing to phase differences between 
the central and outer portions, the rate of 
increase in frequency with gap may also be 
diminished. 

One more possible cause of diminution in slope 
of the characteristic line with increasing gap may 
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Fic. 8. Characteristic lines showing dependence of fre- 
quency upon gap for various crystals. 
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be mentioned. If over the entire gap range the 
primary vibration mode is coupled to a high 
overtone of a secondary mode, this latter mode 
also having a frequency dependent upon the gap, 
then the observed U would be expected to 
diminish with increasing gap. For, in the case of 
the overtone, there are upon each surface of the 
crystal closely adjacent regions of positive and 
negative electrification, and as the gap increases, 
lines of electric intensity tend more and more to 
pass in a lateral direction from one such region 
to its neighbors instead of producing a depolar- 
izing field in the crystal. The effective stiffness 
for this mode, and hence its frequency, is thereby 
decreased. From coupling theory it follows that 
the observed coupled frequency is also slightly 
decreased. That is, the rate of increase of ob- 
served frequency with increasing gap is less than 
if the overtone were absent, and since the effect 
increases with increasing gap, the characteristic 
line must diminish in slope as the gap increases. 

While these considerations suffice to account 
qualitatively for the observed anomalies in the 
gap effect, it remains impossible to assign specific 
reasons for individual cases. It is significant that 
the best agreement with theory, at least for 
small gaps, is found with the BT- and AC-cuts, 
both of which, especially the AC, are inclined at 
such an angle with the z axis as to reduce the 
coupling effects. Only comparatively weak 
evidences of coupling were, in fact, observed 
with these two cuts. 

Of the X-cuts, X 2 had a practically constant 
but very low U at all gaps, over a wider range of 
gap than X 1. Both crystals had high resistances 
and required relatively high voltages for excita- 
tion. At small gaps, X 1 (a circular plate) had 
three prominent response frequencies, separated 
by a few hundred cycles. The results are repre- 
sented by the characteristic lines X 1, X 1’, 
X 1” in Fig. 8. Table III contains only the value 
of U for X 1, for very small and for maximum 
gaps, since this mode was the most pronounced 
of the three. On the whole, these results indicate 
that the almost universal practice of letting 
X-cut plates for thickness vibrations rest hori- 
zontally upon one electrode is not conducive to 
best operation. 

A conclusion that may be drawn from all the 
observations is that a large increase in frequency 
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with gap, approximating the theoretical value, 
is a sign of a good resonator. It is hardly con- 
ceivable that any disturbing circumstance could 
make the observed U any greater than the the- 
oretical value for the primary mode. 


Effect of changing position of crystal in gap 

§20. Brief mention should be made of the 
effect upon frequency of altering the position of 
the crystal between the electrodes. A test was 
made with crystal Q 30 (unsilvered) in which, 
while the total gap was left constant at 1.7 mm, 
the frequency was observed while the bar, always 
parallel to the electrodes, was moved by short 
stages from one side to the other. The frequency 
was a maximum with the bar in the central 
position, decreasing to a minimum as the bar 
was moved toward either electrode, and on still 
closer approach rising to a pronounced maximum, 
followed by a slight decrease just before contact 
with the electrode. The total variation was not 
over 5 cycles per sec., and played a quite insig- 
nificant part in determining the value of U. The 
maximum in frequency for the central position 
was found by Dye (reference 2, p. 450), and also 
still earlier by W. F. Powers in this laboratory. 
Both observers, however, recorded only a 
decrease in frequency as the crystal approached 
either electrode. It seems probable that this 
maximum is due to the finite width of the bar, 
by reason of which, as was explained above, 
some of the lines of electric intensity curve back 
around the edges; this effect, and the consequent 
diminution in frequency, may be expected to be 
slightly more pronounced when the crystal is 
close to one electrode. This view is strengthened 


by the fact that Dye’s bar, which showed a. 


much greater variation in frequency with posi- 
tion between the electrodes than Q 30, was 
relatively less than half as broad. 

The rise in frequency of Q 30 to secondary 
maxima near the electrodes may have been an 
anomalous effect due to some peculiarity in the 
mounting. The voltage was too low to produce 
ionization. As the space between crystal and 
electrode becomes quite small, the periodic com- 
pression of the thin layer of air must increase 
the effective stiffness to a slight extent, but a 
simple calculation shows the effect to amount to 
a fraction of a cycle per second, and in fact the 


secondary maxima did not manifest themselves 
in the experiments in which the total gap was 
varied. Likewise negligible must be the increase 
in frequency due to the slight cooling of the 
vibrating quartz as it approaches an electrode. 
On the other hand, the rise in frequency may 
conceivably have been due to the approach to 
the first air-resonance point on the wide side of 
the gap. The decrease in frequency when the 
crystal was brought extremely close to one elec- 
trode may have been caused by increased air 
friction, as was reported by Watanabe (reference 
3, second citation, p. 710). 

These results point to the great importance, 
in work of high precision, of keeping the crystal 
in a fixed position with respect to the electrodes, 
and of mounting it in a vacuum. 


Observations on Rochelle salt 


§21. For comparison with the gap effect in 
quartz, a word may be added concerning a test 
made by the writer with Rochelle salt. The length- 
wise vibration frequency at zero gap of a Y-cut 
bar 2.5 mm thick, measuring 21.35.9 mm in 
directions at 45° to the x and z axes, was roughly 


56,800 cycles per sec. With a gap of 3.2 mm the 


frequency was 58,000, or about 2 percent greater. 
The effective value of Young’s modulus was 
co=1.03(10") dyne/cm*. On applying Eq. (46) 
and calling the dielectric constant k=10,” we 
find for the piezoelectric constant 62; the value 
134 (10-*). This is in good agreement with the 
value 138 (10-°) recorded by Valasek,** and 
somewhat less than Pockels’ value*® of 165 (10-°). 
These observations indicate that if the elasticity 
and dielectric constant of a piezoelectric crystal 
are known, a rough approximation, at least, for 
the piezoelectric constant can be derived from 
gap observations. 


Results of other observers 


§22. Fig. 17 of Dye’s paper? shows the de- 
pendence of frequency upon gap for thickness 
vibrations of an X-cut quartz plate 15.64 16.40 
X6.15 mm. The writer has computed U from 
Dye’s curve and found a mean value of 0.00536, 
about one-third of the theoretical value. Better 


37H. Miiller, Phys. Rev. 47, 175 (1935). 


38 J]. Valasek, Science 65, 235 (1927). ; 
39 F, Pockels, Abh. d. kgl. Ges. d. Wiss. zu Gottingen 39, 
183 (1894). 
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agreement with theory is found from Fig. 30 in 
Vigoureux (reference 18, p. 68), for an X-cut 
quartz disk 9.43 mm diam., 0.58 mm thick. The 
theoretical value of U for this crystal is 0.0150, 
the value from the curve is 0.0137. Other records 
of frequency-gap observations have been pub- 
lished, for example, by Namba and Matsumura,” 
Koga*® and Watanabe,’ but the data are either 
insufficient or taken with the crystal in the 
generating circuit. Results obtained by the latter 
method are not characteristic of the resonator 
alone. 


VI. OVERTONE VIBRATIONS 
Gap effect for overtones 


§23. Since resonators have sometimes been 
employed with the crystal vibrating at an over- 
tone frequency, it is of interest to examine the 
effect of varying gap on such frequencies. It is 
well known that, if f; is the fundamental fre- 


1. Koga, Proc. I. R. E. J8, 1941 (1930). 
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quency of thickness or lengthwise vibrations, a 
crystal will resonate at a frequency f,=nf; 
approximately, where m is an odd integer. The 
writer has derived expressions for the equivalent 
constants for overtone vibrations, from which 
equations for the gap effect, analogous to Eq. 
(44), have been derived. Only the final results 
can be presented here. 

With thickness vibrations it is found that 
Af n/fn=(1/n*)(4fi/fi), or in other words, the 
relative effect of gap on frequency is only 1/n? 
as great for the overtone as for the fundamental. 
With lengthwise vibrations the value of Af,/f, 
for overtones turns out to be the same as for the 
fundamental. By using an even number of pairs 
of electrodes distributed uniformly along the 
length of the bar, even harmonics (approximately) 
of lengthwise vibrations can also be excited, in 
which case also the gap effect is the same as for 
the fundamental. These relations have not yet 
been tested experimentally. 
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Further Studies on the Mathematical Physics of Metabolizing Systems with 
Reference to Living Cells. II 


N. RASHEvVsky, University of Chicago 
(Received May 6, 1936) 


A system of nonlinear differential equations, established in a previous paper, and describing 
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the diffusion in a metabolizing system, which glycolizes and oxidizes sugar, is approximately 
integrated. With such values of the constants as are suggested by considerations of some 
biological systems, the average approximation obtained by this method is about three to five 
percent. The results are discussed and it is shown that they throw light on the variation of 


oxygen consumption of cells with varying oxygen concentration. 


N connection with a series of previous studies 

of the physics of the metabolizing systems,':?:* 
we have discussed in a preceding paper‘ the 
system of differential equations which describes 
the consumption of sugar with subsequent 
glycolysis and oxidation. Using the same nota- 
tions as before,* the system of equations is as 
follows : 


Vex’ =0, (c) 


with the boundary conditions,‘ for r=ro: 


Dy(de,/dr) (a) 
(b) 

D2(de2/dr) =he(c2’—c2), (2) 
Dz! (des'/dr) =D3(de2/dr), (d) 

D;(dc3/dr) =h3(c3’—cs), (e) 
=D3(de3/dr).  (f) 


We have discussed‘ the solution of this system 
for a very special limiting case, namely that of 
very large D’s. We also have considered for that 
case the possible variations of the permeabilities 
h with the rates of metabolism. In the present 
paper we shall give a more general approximate 
solution of (1) and (2) for the case of spherical 
symmetry, considering however all h’s and D’s 
as constants. 

Eqs. (ia) are independent on the others. For 
spherical symmetry, Eqs. (1a) reduce to 


(d?e1/dr*) +(2/r)(dei/dr) (a) (3) 


d’c,' (b) 


1 N. Rashevsky, Physics 5, 374 (1934). 
2 N. Rashevsky, Physics 6, 35 (1935). 

N. Rashevsky, Physics 6, 117 (1935). 
4N. Rashevsky, Physics 6, 343 (1935). 
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A substitution 
and (4) 
transforms (3) into 
d*u;/dr?=au;/D,; /dr?=0. (5) 


Since a>0*, therefore the general integral of the 
first of Eq. (5) is 


u;=A, sinh (a/D;)'r +B, cosh (a/D,)'r, (6) 


where A, and B, are constants of integration, 
and therefore 


{sinh (a/D;)'r]/r. 


The requirement that c; remain finite for r=0 
makes B,=0. Hence 


c:=A[sinh (a/D,)!r]/r. (7) 


The general integral of the second equation of 
(5) is 


or =A (8) 


Since for ¢:’=Co1, Co. denoting the ex- 
ternal concentration at infinite distance, we have 


(9) 


The remaining two constants A; and By,’ are 
determined by introducing (7) and (9) into (2a) 
and (2b). In this way after some elementary 
calculations the exact solution of (1a) and (2a) is 
obtained, which however is very cumbersome. 
A very good approximation leads however to a 
much simpler expression. We have, namely : 


[sinh (a/D;)'r 
(10) 


With the range of values of a, r and D, which 
we have considered hitherto, and which are 
taken from biological considerations,‘ we find 
that the first two terms of expansion (10) 
constitute a very good approximation. Omission 
of the higher terms introduces an error of only 
one to five percent. We have therefore with a 
good approximation 


(11) 

Hence dc,/dr=ar/3D, (12) 
while (9) gives 

dey’ /dr= —B,'/r’. (13) 


Substituting ro for r in (12) and (13) and intro- 
ducing them into the (2b), we easily find 


Introducing (14) into (9) and the latter 
together with (11) into (2a), we find after simple 
rearrangements 


(15) 


(9), (11), (14) and (15) give us the desired 
approximate solution of (la). _ 

If we introduce (11) into (3a) we shall find 
that the left-hand side or D,V?c,; becomes equal 
to A,a. Since however D,V*c; is equal at any 
point to the true rate of consumption of c; at 
that point, the introduction of the approximation 
(11) is physically equivalent to introducing, in- 
stead of the rate of reaction ac; per unit volume, 
which varies from point to point with ¢, an 
average constant rate of reaction g:=A,a. We 
shall use this to solve approximately the other 
four Eqs. (1). 

The first equations of (1b) and (1c), trans- 
formed to polar coordinates, give 


|= —acit+feecs, (16) 
D;{d?c3/dr?+ (2/r)(de3/dr) |=1.07feocs. (17) 


Solving (17) for fcscz; and introducing it into 
(16), we have: 


(2/r) (dc2/dr) | 
= —ac,+0.94D;[ d2c3/dr* 


+(2/r)(des/dr)]. (18) 
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Introducing for c, the expression (11) and making 
the substitution : 


Co=U2/P; (19) 
we find 
=0.94D3(d?u3/dr’) 
-A /6D, —A (20) 
Integrating this twice, we find 
= 0.94D3u3 —A 
(21) 


a and £6 being constants of integration. 
Introducing again cz and ¢3 according to (19), 
we have 


Doce = 0.94D 3c; /6 
(22) 


Since both c2 and c; must remain finite for r=0, 
we must have 8=0. Therefore (22) gives 


0.94(D3/De)cs3 —A ar? 
—A (23) 


If we knew c3, we would know C2, a being de- 
termined by the boundary conditions. c; is de- 
termined in its turn by (1c). Using the above 
outlined idea, we shall introduce instead of the 
variable rate of reaction for cs, 1.07fcec3, an 
average constant rate G3=1.07fC2c3, defined by 


93=[1/(4/3) are] X1.07f40 f (24) 


If the rate of consumption of c; were constant 
and equal to @;3, then the solution of (1c) with 
(2e) and (2f) would be! 


(25) 


C3’ =Cos— (Garo®/ 3D3’) (1/r), (26) 


Cos being the external concentration of oxygen at 
infinite distance from the system. Hence these 
equations will express approximately the dis- 
tribution of c; and c;’, if we can express Q3 in 
terms of the other constants of the system (1) 
and (2). 

To this end we introduce the abbreviations: 


mM = Co3 — Gst'o/3hs— Gat? ; 


(27) 
n =9;/6D3. 
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(25) thus becomes 
C3=m+nr’. (28) 
Introducing this into (23), we find 


C2= 0.94(D;/Dz) (m +nr*) ar’ 
—A (29) 


Using again the previous values of the various 
constants,‘ we find that the term A ,a’r*/120D,D, 
is small as compared with Aj,ar?/6D. and we 
shall therefore omit it. The maximum error 
introduced by this is less than five percent. 
The average error is about two percent. The 
general integral of the second equation (1b) is 
again of the form:! 


C2’ (30) 


Differentiating (29) and (30) with respect to r, 
substituting ro for r and introducing into (2d), 
we find after simple rearrangements 


By! =(A,a—5.6D3n)r,3/3D (31) 


Introducing this into (30), substituting then ro 
for r in (30) and in (29), and introducing both 
into (2c), we find after some rearrangements 


a= Decor —0.94D 3m (0.94D3n —A 
— (5.6D3n — Aa) De’ (32) 


We have thus expressed a in terms of m and n, 
or, on account of (27), in terms of the yet un- 
known quantity 

Substituting now (27) and (32) into (29), and 
omitting A ,a*r*/120D,Dz, we find after elaborate 
but quite elementary rearrangements 


tA 
+A 
(33) 


On the other hand (25) can be written 


€3=Co3— 


Or, introducing the notations 


+0.16(ro?— 77) / Dae, (35) 
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6= (2D3D3'ro+ 
h3D3;'r*) /6D3D3'hz, 


we have 


C3=Co3— (36) 
Hence, 


— (6co2+ (37) 


The quantities y, 7 and 6 are quadratic functions 
of r. Hence the average 1.07fcsc3 is given by: 


1.07 = 1.07f[co2cos+ 
— (6co2+ (38) 


with 7 =[1/(4/3)ar (39) 


and similar expressions for 7 and 6. 
From (35) 7, 7 and 6 are very easily calculated, 
according to (39), and we find 


/3h2D2' +A are /15D2, 
/h2D2' 
6= [2D3D3'ro+2Dshsre 

|/6D3D3'h3. 


But, according to (24), 1.07fcoc3 is itself equal 
to g3. Hence we have now an equation for the 
determination of g;, by making the right-hand 
side of (38) equal to gs 


(40) 


q3= 1.07f[co2¢03 — 
— ]. (41) 


Eq. (40) expresses 7, 7 and 6 in terms of the 
parameters of the differential equations. Using 
previous values for a, f, 70, /t1, hs and taking’: 
D~10-7, we find that 7 and 6 are of the order 
of 1, while 7 is of the order of 10-*. On the 
other hand @; is nothing else but the average 
rate of oxygen consumption, which is known for 
biological systems considered here to be of the 
order of magnitude of 10-* g-cm~*-min. With 
those values we find that the quadratic terms 
in (41) can be neglected and we find with good 
approximation 


+é6y+7¢03) (42) 


The same result can of course be obtained by 
solving the quadratic equation (41) for g;, and 
then omitting small terms. 
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We now have a complete approximate solution 
of Eqs. (1) and (2). c: and c,’ are completely 
described by (11), (8), (14) and (15) in terms of 
a, f, hi, he, hz, D,, Dy, Ds, D,’, D3, and To. 
In terms of the same parameters c; is described 
by (25), (42) and (40); while (23) and (32) 
give us an expression for Cz. cs’ and c;’ are 
given by (30), (31) and (26), (42), (40) corre- 
spondingly. With the approximation here used 
all internal concentrations are of the form A+ Br’ 
and all external of the form E+G/r, the coeffi- 
cients A, B, E and C being expressed in terms of 
the 12 parameters a, f, hy, etc. 

A measure of the approximation thus obtained 
is given by the amounts by which y, 7 and 6 
‘deviate from their corresponding averages y, 7 
and §. With the above discussed values for the 
physical constants used, a numerical estimation 
shows that the deviation of the extreme values 
from the averages is below three percent. The 
average error is therefore about one to two 
percent. 
Because of a close interrelation between the 


three reactions described by (1), the concentra-- 


tion of any substance depends not only on the 
diffusion coefficient and permeability for this 
substance, but also on the diffusion coefficients 
and permeabilities of the other substances. A dis- 
cussion of those rather complex relations from 
the point of view of biological applications will 
be given elsewhere. Here we shall discuss only 
one general property of our solution, which is 
also of great biological interest. We remark that 
the total rate of oxygen consumption Q; in the 
system is equal by definition to 


Qs= (4/3) (43) 


Eq. (42) shows that if all other parameters are 
kept constant, but the external concentration 
of oxygen, Cos, is varied, G3 is zero for co3=0, and 
with increasing co; increases to an asymptotic 
value 

(co.—)/n. (44) 


Hence Q; will also increase from zero to an 
asymptotic constant value. 
It has been observed® that in living cells the 


°R. W. Gerard and P. S. Tang, J. Cell. Comp. Phys. 1, 
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total oxygen consumption per cell is constant 
for high oxygen concentrations and drops to 
zero for low oxygen concentrations. An interpre- 
tation of this phenomenon has been attempted 
by R. W. Gerard® on the assumption of a strictly 
constant true rate of oxygen consumption, and 
ascribed to a limitation of the depth of penetra- 
tion of oxygen due to a finite diffusion coefficient. 
A more general treatment has been given by N. 
Rashevsky,’ so as to include effects of perme- 
ability. The assumption of a constant rate of 
reaction, independent of concentration, is how- 
ever surely untenable for biological systems. 
The present study, free of this limitation, gives 
us a much more general insight into this phe- 
nomenon. In the present case it is due largely to 
the circumstance that the true rate of oxygen 
consumption is proportional to Csc3, that is, to 
both the concentration of oxygen and that of 
lactic acid. An increase of oxygen concentration 
results in an increased consumption both of 
oxygen and of lactic acid. This decreases Co. 
The product c2c3 increases therefore less rapidly 
than c; alone, and tends to a constant value for 
infinite c3, C2 tending to zero. 

This effect is found also when using the 
quadratic equation (41), which involves more 
calculations. The asymptotic value (44) for @; is 
still in the range for which the quadratic term 
may be neglected. 

Eq. (31) shows the interesting effect of the 
glycolytic constant a. If a is large enough, By’ is 
positive, the concentration of lactic acid outside 
the system is decreasing with increasing 7, 
which corresponds to a net production of lactic 
acid. This is the usual case. If a is very small, 
B,’ becomes negative. Lactic acid flows into the 
system. This could happen, when the breakdown 
of sugar is not rapid enough to supply sufficient 
lactic acid for oxidation. In the first case co2 may 
be made zero, without impairing the whole 
system of reactions. In the second case ¢o2 must 
be kept finite. In order for the reactions to take 
place at all, there must be a supply of lactic acid. 


503 (1932); W. R. Amberson, Biol. Bull. 55, 79 (1928); 
C. S. Shoup, J. Gen. Physiol. 13, 27 (1929). 
®R. W. Gerard, Biol. Bull. 60, 245 (1931). 
7 N. Rashevsky, Protoplasma 20, 125 (1933). 
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